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Abstract 

We investigate the thermal equilibrium properties of kinks in a classical 4 field 
theory in 1 + 1 dimensions. The distribution function, kink density, and correlation 
function are determined from large scale simulations. A dilute gas description of 
kinks is shown to be valid below a characteristic temperature. A double Gaussian 
approximation to evaluate the eigenvalues of the transfer operator enables us to extend 
the theoretical analysis to higher temperatures where the dilute gas approximation 
fails. This approach accurately predicts the temperature at which the kink description 
breaks down. 
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The statistical mechanics of coherent structures such as solitons and solitary waves 
in nonlinear systems has been a subject of study for some time |l|. Recent interest has 
been fueled by new applications not only in condensed matter physics [@||, but also 
by potential applications in particle physics (sphalerons) |4j] and cosmology (domain 
walls, baryogenesis) f5|. In this paper we focus on the classical equilibrium statistical 
mechanics of solitary wave solutions ( "kinks" ) of a tachyonic mass </> 4 field theory in 
1 + 1 spacetime dimensions with Lagrangian 

L = \ - \ (d x ®) 2 + \m 2 & - \\^. (1) 

For our simulations we use the dimensionless form of this theory, given by the trans- 
formations: $ — > acj), x — > x/m, and t — > t/m, where a 2 = m 2 /\. The equation of 
motion then becomes 

d 2 tt <t> = d 2 x J - <P (0 2 - l) . (2) 

What makes this model so useful is that its behavior is representative of a large class 
of soliton-bearing systems. Moreover, it is amenable to both theoretical analysis and 
numerical simulation. 

The statistical mechanics of kinks in this system has been studied by two ap- 
proaches. In the first, and phenomenological, approach one assumes that the kinks 
and the fluctuations ( "phonons" ) about the asymptotic field minima may be treated 
as weakly interacting elementary excitations. The canonical partition function can 
then be found by standard methods JT| |J . Alternatively, as shown by Krumhansl 
and Schrieffer (KS) ||, building on earlier work of Sears, Scalapino, and Ferrell [|J, 
it is possible to calculate the partition function, in principle exactly, by exploiting a 
transfer operator technique. KS showed that in the low temperature ("dilute gas") 
limit the partition function factorizes into a contribution from a harmonic term and 
from a tunneling term which they were able to identify with the phonon and kink 
contributions respectively in the phenomenological theory. The ideas of KS were fur- 
ther refined and extended to a wider class of systems by Currie, Krumhansl, Bishop, 
and Trullinger 0. In particular, interactions of kinks with linearized phonons were 
considered, leading to substantial corrections to the results of KS. 

A key result of these efforts is the prediction that the spatial density of kinks 

n k oc yfE k f3exp(-E k f3), (3) 

where E k = ^8 /9m 3 / A (E k = ^8/9 for the dimensionless form of the theory) is 
the kink energy, and j3, the inverse temperature (for the dimensionless case, (3 — > 
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(3 j (a 3 vA)) • A related prediction is that at low temperatures the field correlation 
length A defined by 

(#OMaO)~exp(-a?/A) (4) 
has an exponential temperature dependence [[/J, 

Computer simulations to verify these results date back to Koehler, Bishop, Krumhansl, 



and Schrieffer Q who found only a qualitative agreement. Recent work |L0| |TT| [12| has 
led to more detailed comparisons, however significant discrepancies have been found 
over the temperature range explored in these simulations. It has been speculated that 
these discrepancies are due to finite size effects and phonon dressing of the bare kink 
energy (breather contributions to the free energy may also be significant ||13|| ). In this 
paper we show that the difficulty is partly that previous simulations were not carried 
out at low enough temperatures. Indeed, for the range of temperatures studied in 
Refs. [|K]] QTTJ |R] , the dilute gas approach is simply not a valid approximation (there 
is also a problem with the operational definition of kinks used by these authors over 
the range of temperatures they studied). 

We have studied the equilibrium statistical mechanics of kinks in the 4 model (|I]) 
by implementing a Langevin code on a massively parallel computer. The key idea is 
to supplement the equation of motion (0) with noise and viscosity terms obeying an 
appropriate fluctuation-dissipation theorem so that the system is driven to thermal 
equilibrium at the desired temperature |14|| . To understand our results in the high and 
intermediate temperature region not susceptible to a dilute gas analysis, we have used 
a double Gaussian wave function approximation for the quantum mechanical problem 
which results from applying the transfer operator method. Our results are: (1) the 
dilute gas predictions for the kink density and the correlation length are very accurate 
below a certain (theoretically estimable) temperature, (2) above this temperature the 
double Gaussian results for the kink number and the correlation length agree with the 
simulations, (3) kinks are found to disappear above a characteristic temperature, in 
good agreement with our theoretical prediction, and (4) our Gaussian approximation 
accurately describes the classical single point field distribution function at high and 
intermediate temperatures where the dilute gas (WKB) approximation breaks down. 
We also observe a temperature dependence of the asymptotic field value away from 
a kink which may be attributed to kink-kink interactions. 

The canonical partition function for the Lagrangian ([!]) is given by the functional 
integral 

Z= f D<f>Dir exp 7r)) (6) 



where 7r is the canonical momentum of the field and H, the field Hamiltonian. The 
transfer operator technique |§ reduces the calculation of the partition function in the 
thermodynamic limit to simply finding the ground state energy of the double well 
quantum Hamiltonian (here written for the dimensionless case), 

where = (5<f>. At low temperatures the two wells are widely separated and the 
ground state energy is given by the oscillator ground state energy for one of the wells 
minus the tunnel-splitting term, usually calculated by WKB methods. The dilute 
gas/WKB approximation for the kink number is valid when the tunnel-splitting is 
small enough such that only the first two energy eigenstates are necessary to estimate 
the ground state energy of the Hamiltonian (|7|) ||15|| . At higher temperatures where 
kinks still exist, higher energy states cannot be ignored. Since kinks are associated 
with tunneling, we expect them to vanish when the ground state energy is higher 
than the classical barrier height: this intuition is confirmed by our simulations. 

One can compare the simulations of the kink system with numerical solutions for 
the energy eigenvalues of the Hamiltonian Hq. Instead, we take a different approach 
by implementing a double Gaussian variational method (see [TBI for details) which is 
an order of magnitude more accurate than the simple Gaussian approximation ]l7j 
for this problem and correctly accounts for the reduction of energy due to overlap 
terms in the wave functions of the two wells, at least for moderate to large overlaps. 
Three qualitatively different regimes exist: (1) all the energy eigenvalues lie above the 
classical barrier, (2) the ground state energy lies below the classical barrier height, and 
(3) the energy difference between the ground and first excited state becomes negligible 
in comparison with the energy difference between the ground and the second excited 
state (this occurs at /3 ~ 6.7). Our simulations confirm the theoretical expectations 
that there are no kinks in region (1), that there are kinks, but that the dilute gas 
approximation is invalid in region (2), and finally, that the dilute gas approximation 
is accurate in region (3) (a regime unexplored in detail by previous simulations). 

We measure the classical single point field distribution function P[<p] directly in 
our simulations. For the analogous quantum mechanical problem this is just the 
square of the ground state wave function ^/ - Results from the simulations and our 
theory are compared in Fig. I and are in reasonable agreement. The presence of kinks 
implies a double peak in P[4>] ||18||(the converse is false) while a single peak at the 
origin means that kinks and thermal phonons can no longer be distinguished. From 
the simulations such a transition occurs at /3 1.7, in agreement with the theoretical 
calculation of when \l/ 2 goes over from a double to single peaked distribution. As 
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expected, this is also the temperature (/3 = 1.734) where the ground state energy 
crosses the classical barrier height (a discussion of various methods to determine the 
characteristic temperature is given in Ref. [|19|). The double peaks in the distribution 
function move inward from the classical minimum as the temperature increases (this 
is clearly seen in Fig. I). We attribute this to kink- kink interactions with the kinks 
presumably sacrificing some potential energy to lower their strain energy. 
The Langevin equation for the dimensionless theory is 



To guarantee an approach to equilibrium, the Gaussian, white noise F and the vis- 
cosity 7] are related via the fluctuation-dissipation theorem: 



We carried out numerical simulations on lattices with 16384 sites and solved the 
Langevin equation using a standard leapfrog algorithm. The time step was At = 0.02 
and the lattice spacing was Ax = 0.5. For all of the simulation results reported here 
we used a viscosity of r\ = 1.0. 

Our system size is one to two orders of magnitude larger than that in most previ- 
ous simulations. Large system sizes are necessary to get acceptable statistics at low 
temperatures. Systems were evolved from a random initial condition to equilibrium. 
The length of time necessary to ensure equilibrium increased with inverse temper- 
ature. For (3 = 8 the time required was approximately 10 7 time steps, and for the 
highest temperatures, less than 10 5 steps. 

Two quantities of interest reported here are the kink number and the field corre- 
lation length. To compute the kink number, we need an operational way to identify 
kinks, although there is an exact kink solution available theoretically. As a result we 
examine several possible definitions, all of which rely on a knowledge of the kink size. 
From the classical solution for a kink centered at xo, 4> = tanh((x — xq)/ v^), the kink 
scale Lk is approximately 8 lattice units. Raw kink configurations are shown in Fig. 
II. At low temperatures (/3 > 6), kinks may be identified easily, however at higher 
temperatures this is clearly not the case. 

The simplest thing to do is to count the number of zero-crossings of the field, 
since one may argue that these are the "tunneling events" which correspond to kinks. 
However, at higher temperatures there are zero-crossings due to thermal noise, and 
counting all zero-crossings would lead to a gross overestimation of the number of 
kinks. A possible solution is to use a smoothed field by either averaging or block- 
spinning the actual field configuration over a length of the order of the kink scale. 






(9) 
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The latter approach was taken in previous simulations [1C][5 ]|[12 . This solution is 
not without flaws either, as rapid fluctuations can still appear as kinks. We prefer to 
count kinks in the following way: at a particular time we first find all zero-crossings. 
To test the legitimacy of a given zero-crossing we check for zero-crossings one kink 
scale (8 lattice units) to its right and to its left. If no zero-crossings are found, we 
count it as a kink, otherwise not. 

The number of kinks is plotted against (3 in Fig. III. Above (3 ~ 6, the aver- 
aged field method and our method for counting kinks agree. Moreover, in this (low 
temperature) range, the dilute gas expression for the kink number ([|) is in excellent 
agreement with the data. At elevated temperatures, there is a clear disagreement 
between the two methods of counting kinks. The average field technique has the 
number of kinks monotonically increasing with temperature; whereas, in accord with 
intuition and the behavior of P[4>], the second technique clearly shows a reduction 
in the kink number at higher temperatures. Moreover, in this temperature regime 
the number of kinks computed with the smoothing method depends strongly on the 
smoothing scale. We conclude that for (3 < 6, the number of kinks cannot be extracted 
with any confidence from the smoothing method. Unfortunately, this is precisely the 
temperature regime explored in previous simulations. 

The correlation length A is plotted against (3 in Fig. IV. For (3 > 6, the WKB 
prediction (|J) holds, whereas for (3 < 5, where the wave function overlaps are not 
negligible, the double Gaussian approximation is valid. Fortunately, there are no 
ambiguities here with regard to measurements at higher temperatures. 

As a consequence of the above results, we conclude that the dilute gas/WKB 
approximation is excellent for (3 > 6 with no further phonon dressing of the bare kink 
energy beyond that already included in (||) and (||) at these low temperatures. At 
higher temperatures, the WKB analysis fails, though theoretical progress is possible 
with the double Gaussian technique. 
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Figure Captions 

Fig. I: The classical distribution function P[<f>] given by the simulation (solid line) and 
the distribution \l/ 2 from the double Gaussian approximation (dashed line) plotted 
against <j> for (3 = 2. The potential V = -(1 - <j) 2 /{2f3 2 ))<f) 2 /{2(3 2 ) of the equivalent 
quantum problem (0) is also shown. 

Figure II: Field configurations, from top to bottom, at (3 = 2, (3 = 4, and (3 = 8. 
Only a 1000 lattice unit sample of the total lattice size of 16384 is shown. 
Figure III: Total number of kinks and anti-kinks as a function of (3. Squares denote 
counts with a smoothed field (smoothing scale=8 lattice units) definition of kinks, 
triangles for the zero-crossing count method explained in the text, and the solid line 
is a fit to (Q). Where not shown, error bars are of the order of the sizes of the symbols. 
The smoothed field result is in agreement with the results of Ref. |T2[ for (3 < 5. 
Figure IV: Field correlation length A as a function of (3. The short-dashed line is the 
Gaussian prediction, whereas the long-dashed line is the WKB result (|5|). 
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